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Abstract 
O 

u 

£^j We develop an optimized perturbation theory for the Ginzburg - Landau 

description of thermal fluctuations effects in the vortex liquids. Unlike the 
high temperature expansion which is asymptotic, the optimized expansion is 
convergent. Radius of convergence on the lowest Landau level is ax = — 3 
in 2D and ay = —5 in 3D. It allows a systematic calculation of magnetiza- 

^h \ tion and specific heat contributions due to thermal fluctuations of vortices 

in strongly type II superconductors to a very high precision. The results are 
in good agreement with existing Monte Carlo simulations and experiments. 
Limitations of various nonperturbative and phenomenological approaches are 
noted. In particular we show that there is no exact intersection point of the 
magnetization curves both in 2D and 3D. 
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I. INTRODUCTION 

Thermal fluctuations play a much larger role in high T c superconductors than in the low 
temperature ones because the Ginzburg parameter Gi characterizing fluctuations is much 
larger |l[ . In addition the presence of magnetic field and strong anisotropy in 
superconductors like BSCCO effectively reduces their dimensionality thereby further 
enhancing effects of thermal fluctuations. Under these circumstances the mean field line 
separating Abrikosov lattice from "normal" phase becomes a phase transition between 
vortex lattice and liquid far below the mean field phase transition line [@,|l| clearly seen in 
both magnetization [Q] and specific heat experiments [|J . Between the mean field transition 
line and the melting point physical quantities like the magnetization, conductivity and 
specific heat depend strongly on fluctuations. Several experimental observations call for a 
refined precise theory. For example a striking feature of magnetization curves intersecting 
at the same point (T*, H*) was observed in a wide rage of magnetic fields in both layered 
(2D or quasi 2D) materials and more isotropic ones 0. To develop a quantitative 
theory of these fluctuations even in the case of the lowest Landau level (LLL) 
corresponding to regions of the phase diagram "close" to H C 2, is a very nontrivial task and 
several different approaches were developed. 

Long time ago Thouless and Ruggeri [|7]|| proposed a perturbative expansion around a 
homogeneous (liquid) state in which all the "bubble" diagrams (see Fig. 5) are resummed. 
Unfortunately they proved that the series are asymptotic and although first few terms 
provide accurate results at very high temperatures, the series become inapplicable for LLL 
dimensionless temperature ar ~ (T — T m f(H))/{TH) 1 / 2 smaller than 2 in 2D quite far 
above the melting line (believed to be located around ax = — 12). Generally attempts to 
extend the theory to lower temperatures by the Borel transform or Pade extrapolation 
were not successful ||. Several nonperturbative methods have been also attempted. 
Originally the RG method was proposed || and developed |I| although, since the 
transition is first order, no solutions of the RG equations can been found. The set of 



perturbative "parquet" diagrams [II J have been resummed and the large N limit have 



been considered [|12|. Tesanovic and coworkers developed a method based on an 



approximate separation of the two energy scales Jl3| in both 2D and 3D. The larger 



contribution (98%) is the condensation energy, while the smaller one (2%) describes 
motion of the vortices. The theory explains the intersection of the magnetization curves. 



This question has been tackled in 2D by rather phenomenological approach in [ll] . Some 
Monte Carlo simulations are available |H)| , |T6|| . 

Meantime experimental precision increased dramatically. New methods like measurement 
of magnetization using the Hall probes || were invented. One can achieve a precision that 
allows clearly to see a tiny magnetization jumps of only O.lOe in BSCCO and a sharp peak 
in specific heat in YBCO. 

In this paper we apply optimized perturbation theory (OPT) first developed in field theory 
|H| |T9"| to both the 2D and 3D LLL model. It allows to obtain a convergent series (rather 
than asymptotic) and therefore to calculate magnetization and specific heat of vortex 
liquids with definite precision. The precision for various values of the LLL scaled 
temperature a^ are given in Table 3 and 4. The radius of convergence is ay = —3 in 2D 
and a? = —5 in 3D. One the basis of this one can make several definitive qualitative 



conclusions. The intersection of the magnetization lines in only approximate not only in 
3D (the result already observed in Monte Carlo simulation [f[ql ), but also in 2D. The 



theory by Tesanovic et al [HJ in 2D describes the physics remarkably well in high 
temperatures, but deviates on the 5-10% precision level at a T = —2. 

The paper is organized as follows. The models are defined and the general OPT described 
in section II. The 2D and the 3D calculations are described in III. Results and comparison 
with other theories and experiments are given in section IV. We conclude in section V. 



II. MODELS AND GENERAL IDEA OF OPTIMIZED PERTURBATION 

THEORY. 

A. The 2D model 



Our starting point is the Ginzburg-Landau free energy: 
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(By, 0) describes a nonfluctuating constant magnetic field in Landau gauge and 
2 " A tU ~ '" For strongly type II superconductors like the high T c cuprates 
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where A 
D 

(k ~ 100) and not too far from H c2 (this is the range of interest in this paper, for the 

detailed discussion of the range of applicability see |2"0fl ) magnetic field is homogeneous to a 

high degree due to superposition from many vortices. For simplicity we assume 

a(T) = aT c (l — t), t = T/T c , although this temperature dependence can be easily modified 

to better describe the experimental H c2 (T). 

Throughout most of the paper will use the coherence length £ = \Jh 2 / (2m a i,aT c ) as a unit 

dH c2 (T c ) rp 



of length and — -jj, — j. c 
is rescaled as i\) 2 — > %r^ 



2%2 as a unit of magnetic field. After the order parameter field 



ip , the dimensionless free energy (the Boltzmann factor) is: 
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where L z is width of the sample. The dimensionless coefficient describing the strength of 
fluctuations is 



r^TT- 2 + rn ab b' 1 / 327TgVg7; 



(3) 



where Gi is the Ginzburg number in 2D . When 



l-t-b 
12b 



<< 1, the lowest Landau level 



(LLL) approximation can be used [2D]. The model then simplifies due to the LLL 



constraint, —R^ip = \^ to: 
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This reduced model exhibits the LLL scaling. Rescaling again x — > x/yb, 



y^y 



jyfb and 



I 2 -* |-0| 2 -v/ S' one obtains 



J 4vr 



where the 2D LLL reduced temperature 
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buo 2 
is the only parameter in the theory |2l|,[7||. In total, we have done the rescalin^ 



(6) 
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x -> ^x/a/6, y -> ^//v 7 ^ 



(7) 



We will be interested in thermodynamic properties of the model determined by partition 
function Z = j DipDipexp (— /) and will mainly study only the rescaled partition function 
Z r {aT) = f Dtl) r Dijj r exp (— /) = Z/ J, where J is a Jacobian. Consequently to obtain, for 
example, the free energy density from the corresponding quantity in the rescaled model 

feff : 



■ - °f r , one should use the following relation: 
TlogZ T V' (-Air\ogZ r J) 
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From now on we work with rescaled quantities only and related them to measured 
quantities in section IV. 

B. The 3D model 



(8) 



For 3D, the GL model takes a form 

h 2 



d x- 



te 



2rn ab 
which can be again rescaled into 

f = ^ = ~ fd 3 x 



v -Tc A ^ 



+ ^\d^\ 2 + a\^ + b -\^ 
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(9) 



(10) 



by x — > £x, y — > £y, z — > -I72, W ~ *■ -fr^ 2 , where 7 = m c /m a b is anisotropy. The Ginzburg 
number is now given by 



1 /327re^^7_^ 

2 I c 2 /i 2 



(11) 



Within the LLL approximation: 



/ = - = -/ <'"■>■ 
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It also possesses an LLL scaling different from the 2D one. After a rescaling 



f bu) 
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x — > x/yb, y — > y/yb, z — » z ( T~m ) > ^ ~ * ( 4^73 ) ^ 2 ' ^ ne dimensionless free energy- 
becomes: 



2/3 
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AttV2 
The 3D reduced temperature is: 
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The relation between the original and scaled quantity (the 3D Jacobian contains an 
ultraviolet divergent term which will be will cancel with corresponding one loop divergence 
and is not written here) is: 

TlogZ T V (-4nV2logZ r j) 
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C. General description of the optimized gaussian perturbation theory for scalar 

fields. 



We will use a variant of OPT, the optimized gaussian series |19] to study the vortex liquid. 
It is based on the "principle of minimal sensitivity" idea ||17|| , first introduced in quantum 
mechanics. Any perturbation theory starts from dividing the Hamiltonian into a solvable 
"large" part and a perturbation. Since we can solve any quadratic Hamiltonian we have a 
freedom to choose "the best" such quadratic part. Quite generally such an optimization 
converts an asymptotic series into a convergent one (see a comprehensive discussion, 



references and a proof in [IS]). Here we describe the implementation of the OPT idea 
using a simple model of a real scalar fields <ft 

f = -\<PD- l <P + V(<P). 



(16) 



The free energy is divided into the "large" quadratic part and a perturbation introducing 
variational parameter function G^ 1 : 



/ = K + av, K = i^TV, v = f- l -<PG- l <p. 



(17) 



Here the auxiliary parameter a was introduced to generate a perturbation theory. It will 
be set to one at the end of calculation. Expanding the logarithm of the statistical sum to 
order a n+1 



Z= I £>0exp(— K) exp(— av) = j T>4>/ — (av )'exp(— K), 



(18) 



f n {G] = -\ogZ = -\ot 



Vcj)exp(-K) 



i=0 

ra+1 

E 



-a) 



■v. 



('••) 



K ' 



where () K denotes the sum of all the connected Feynman diagrams with G as a propagator 
and then taking a — > 1, we obtain a functional of G. To define the n th order OPT 
approximant f n one minimizes f n [G] with respect to G: 



f n = mm f n [G]. 



(19) 



Till now the method has been applied and comprehensively investigated in quantum 
mechanics only ( pi and references therein) although attempts in field theory have been 
made |L7| . 



III. OPT IN THE GINZBURG - LANDAU MODEL 

A. 2D 

Due to the translational symmetry of the vortex liquid there is only one variational 
parameter, e, in the free energy defined by: 



(20) 



where an = clt — £■ It is convenient to use the quasi - momentum eigenfunctions similar to 
those used extensively in the vortex lattice: 
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where a a == \/^- We expand 



ip(x) 






rf(k) 



Then the propagator in the quasi - momentum basis is: 

mk)m) = —6(k+i) 



(22) 



(23) 



In the coordinate space: 



(i}*(x 1 ,y 1 )ip(x 2 ,y 2 )) 



exp 



Oi - x 2 ) {yi + V2) 



X 



(24) 



exp < -- [(xi - x 2 f + (2/1 - y 2 f] 



The Feynman rules are given on Fig. 1. We have a propagator denoted by a directed line, 
Fig. la, connecting two points (xi,y±) to (x 2 ,y 2 ). For the first term in v, we have a vertex 
represented by a dot on a line, Fig. lc with a value of f-a>H- The second term is a four line 
vertex, Fig. lb, with a value of f-h To calculate the effective energy density 
f e ff = — 47r In Z, we draw all the connected vacuum diagrams. Then one of the coordinates 
is fixed, and all the other are integrated out. We calculated directly diagrams up to the 
three loop order shown on Fig.2,3,4 with the following result. 



h = /o - -^ (18 + 8a H e + a 2 H e 2 ) 

h = h + A ( 662 + 324a ^ £ + 54a2 H £2 + 3a H^ 3 ) 



(25) 



However to take advantage of the existing long series of the non optimized gaussian 
expansion, we found a relation of the OPE to these series. Originally Thouless and 
Ruggeri calculated these series f e ff to sixth order, but it was subsequently extended to 
12 i/l by Hikami et al |22| and to 13 t/l by Hu and MacDonald p3| |. It can be presented using 
variable x introduced by Thouless and Ruggeri j7| 
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n=l 



as follows: 

00 
feff = 21og^ + 2f 2D (x) , f 2D (x) = Y^ c < 

The coefficients are given in Table 1. 

Table 1. 

Coefficient c and z in 2D. 



x . 



(27) 



n 


Cn 


Z n -1 


1 


-2 


-4 


2 


-1 


-6 


3 


38 
9 


-12.239721181139888 


4 


-39- ^ 


-7.508888400035477 


5 


471.39659451659446 


-7.349933383279474 


6 


-6471.5625749551446 


-14.152646217045422 


7 


101279.32784597063 


-9.961364397930787 


8 


-1779798.7875947522 


-9.174960576928443 


9 


34709019.614363678 


-15.232548389083844 


10 


-744093435.66822231 


-11.629924499110746 


11 


17399454123.559521 


-10.8399817525306 


12 


-440863989257.28510 


-15.9366927661989 


13 


12035432945204.531 


-12.753308785106007 



We can obtain all the OPT diagrams which do not appear in the gaussian theory by 
insertions of bubbles and vertex Figlc. insertions from the diagrams contributing to the 
nonoptimized theory. Bubbles or "cacti" diagrams, see Fig.5 are effectively inserted in 
eq.(|2TD by technique known in field theory |24j| : 



feff = 2log^ + 2f 2D (x) 



ft 
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Summing up all the insertions of the mass vertex is achieved by 



62 = s + aajj. 



(28) 



(29) 



We then expand f e ff to order a n+1 , and then taking a = 1, to obtain / n . Calculating / 
that way, we checked that indeed the first three orders agree with the calculation 
performed by a direct calculation. Here a few more terms are displayed: 
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The nth OPT approximant f n is obtained by minimization of f n (e) with respect to e: 

d d 



de da H 



f n (e, a H ) = 0. 



(31) 



The above equation is equal to 2 ^ +3 times a polynomial g n (z) of order n in z = e ■ a H . 
That eq. (|3l|) is of this type can be seen by noting that the function / depends on 



combination 



only. We were unable to prove this rigorously, but have checked it to 



(e+aa H ) 

the 40 t/l order in a. This property simplifies greatly the task: one has to find roots of 
polynomials rather than solving transcendental equations. There are n (real or complex) 
solutions for g n (z) = 0. However (as in the case of anharmonic oscillator |Tj|) the best 
root is the real root with the smallest absolute value,. The roots z n for n = to n = 12 are 
given in Table 1. 

aT + y/af 



We then obtain e(ar) 



t^—Az. 



- solving z n = e ■ an = £clt — £ 2 - For z 



-4,we obtain 



the gaussian result, dashed line marked "TO" on Fig. 6. 



B. 3D 



In the 3D, the LLL Ginzburg - Landau model, we set: 



K 



u 



AnV2 V 

a 
4tt^2 



1 



Off 



+ 2 1 



and 



ip(x) 



k:>, 



exp [izk 3 ] yg k (x) 
(V2T) 3 



1>{k) 



(32) 



(33) 



The propagator is 
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or in the coordinate space: 
(ip(x 1: y 1 z 1 )ip(x 2: y2,z 2 )) 
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Thus the propagator in the coordinate space factorizes into a function of coordinates (x, y) 
perpendicular to magnetic field and a function of the coordinate z parallel to it. The mass 
insertion vertex, Fig. lc, now has a value of TT/ 2 a H, while the four line vertex is tt-^ts- The 
calculation is basically the same as in 2D, the only difference being extra integrations over 
k 3 . However since the propagator factorizes, these integrations can be reduced to 
corresponding integrations in quantum mechanics of the anharmonic oscillator [0.[[^ . 



Again we can take an advantage of existing long series of the nonoptimized gaussian 
expansion [0,^]. The results to seventh order are: 



feff 
hD (x) 



4^ + 4v^/ 3 d (x) ; 
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where v E is given by a solution of the cubic gap equation (y^) — a T ^ — 4 = 0: 

1/3 -, / . n 1/3 



^ = a T [ 54 + 3i/324-3a| 



and coefficient c„ are listed in Table 2: 



+ ^ ( 54 + 3 A /324 - 34 
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Table 2. 

Coefficients c and z in 3D 



n 


Cn 


Z n -1 


f 


-2 


-4 


2 


-0.5 


-5 


3 


1.583333333 


-8.80317864821579 


4 


-12.667361111 


-6.187603657880674 


5 


125.59552619 


-5.960012621607176 


6 


-1430.5928959 


-9.47212746817198 


7 


18342.765997 


-7.430474107869646 


8 


-261118.67703 


-6.907260317913621 


9 


4084812.307 


-9.8195351835546 



Similarly the OPT formula for the effective energy density can be obtained by using the 
generational function: 

a 



feff = 4V / ^1 + 4^/30 (x) , X = A 

2 (vO 
and ^fe{ is given by a solution of equation 

(v^i) 3 - £2v^i - 4a = 
with E2 = £ + actH- The solution of eq.(B9T) can be obtained perturbatively in a: 
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C 2 4 C 2 

Expanding / e // in a to order n + 1, then one then sets a = 1 to obtain /„ 
We list here first few OPT approximants f n 

~ t r 2a H 4 
/o = 4v^ + — £ + - 



e e 



/1 — /o 
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907 + 510a H ^ + 96a 2 H e + 6a 3 H Ve 



(41) 
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228.8335069417501 151.166666666a H 
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The OPT n order result f n (ar) is obtained optimizing f n by varying e 



d_ 
de 



d 



dan 



/„ (e, a H ) = 0. 



(42) 



Similarly to eq.(|31"D in 2D this is equal to 3 i — g n (z) , where now z = clh\/~£ and g n (z) is a 

rank n polynomial. Solving $„ (2) and choosing a real root with the smallest absolute value 
19fl , we obtain z n listed in table 2 up to n = 8. Then we solve for y/e the equation 
= au\[£ = (ot — £ ) a/s- The solution is 



Ve = 2 1/3 a T f -27^ 



\ -1/3 

-108«3 + 7292 2 + 



32V3 



-27^ + 



■108<4 



7292^ 



1/3 



(43) 



IV. RESULTS AND COMPARISON WITH OTHER THEORIES AND 

EXPERIMENTS 

A. Energy, precision of OPT 

On Fig. 6 we present OPT for orders n = (gaussian), 1, 3, 4, 5, 6, 8, 9, 12 together with 
several orders (TO, ...12) of the nonoptimized high temperature expansion in 2D . The 
values of free energy of 2D and 3D models for several a^ are tabulated in Table 3 and 
Table 4 respectively. One clearly observes that in 2D the OPT series converge above 
cit = —2.5 and diverge below a? = —3.5. On the other hand, the non optimized series 
never converge despite the fact that above ay = 2 first few approximants provide a quite 
precise estimate consistent with OPT. Above a^ = 4 the liquid becomes essentially a 
normal metal and fluctuations effects are negligible (see Fig. 7, 8) and are hard to 
measure. Therefore the information the OPT provides is essential to compare with 
experiments on magnetization and specific heat. 

Table 3. 

Free energy f n at different orders (modula constant — 21og47r 2 ) 



or 


-2 


-1.5 


-1 


-0.5 


/o 


-2.19416 


-1.42941 


-0.749027 


-0.146255 


/1 


-2.77516 


-1.80556 


-0.988706 


-0.297222 


h 


-2.53854 


-1.68294 


-0.925643 


-0.264857 


U 


-2.55889 


-1.69143 


-0.92912 


-0.266258 


h 


-2.70076 


-1.74015 


-0.945544 


-0.271734 


h 


-2.62447 


-1.71822 


-0.939384 


-0.270031 


h 


-2.51533 


-1.6923 


-0.933365 


-0.268653 


/10 


-2.59943 


-1.70944 


-0.936772 


-0.269318 


/12 


-2.72613 


-1.73113 


-0.940395 


-0.269915 



11 



If precision is defined as (f 12 - f w ) /f 10 , we obtain 4.87%, 1.27%, 0.387%, 0.222%, 0.032% 

at a T = —2, —1.5, —1, —0.5, respectively. We chose approximants n = 

, 1, 3, 4, 6, 7, 9, 10, 12 because they are "the best roots" in a sense defined in ref. 



19 



section 5. For comparison with other theories and experiments on Fig. 7 and 8 we use the 
10 i/l approximant. 

In 3D the picture is much the same. The series converge above ax = —4.5 and diverge 
below ot = —5.5. The non optimized series are useful only above a? — — 1. 

Table 4. 

Free energy f n at different orders for 3D 



ax 


-5 


-3 


-1.5 


-1 


/o 


-4.73313 





2.65763 


3.41112 


h 


-6.493 


-0.375697 


2.53901 


3.32829 


h 


-6.92585 


-0.427383 


2.5287 


3.3222 


h 


-5.27595 


-0.280923 


2.55551 


3.338 


h 


-5.68059 


-0.292936 


2.55455 


3.33757 


h 


-4.68076 


-0.265834 


2.55647 


3.33839 


h 


-7.32654 


-0.313048 


2.55364 


3.33722 


h 


-5.33149 


-0.301797 


2.55392 


3.33731 


h 


-8.01907 


-0.316175 


2.55359 


3.3372 



We define the precision as (/ 7 — / 4 ) / f 7 . / 4 and f 7 are the best roots among the sequences. 
Then we obtain 6. 55%, 2. 94%, 0.0247%, 0.00779222%, at a T = -5, -3, -1.5, -1 
respectively. 

B. Other theories 

We compare with other theoretical treatments of the same model. A direct method is the 
Monte Carlo simulation of the same model. The 2D model was simulated by Moore, Kato 
and Nagaosa , and Hu McDonald. The circles on Fig. 8 for specific heat are the results of 
the Monte Carlo simulation of the LLL system by Kato and Nagaosa in ref. [15| performed 
with 256 vortices. In 3D the model was simulated with 100 vortices by Sasik and Stroud 
||16|| , magnetization data are compared with our results on Fig. 9. 

An analytic theory used successfully to fit the magnetization and the specific heat data 
p5| was developed in [JI3"]. Their free energy density is: 



feff 

u 



a T U 2 otU 



4 
1 



U 2 a% 



2 + 2arc sinh 



arU 



(44) 



+ 



V2 VP- 



+ tanh 



a T 1 

i72 + 2 



1 
72 



% 



The corresponding magnetization and specific heat are shown as a dashed lines on Fig. 7 
and 8 respectively. The theory applies not only to the liquid phase, but also to the solid 
although the transition is not seen (should be considered as a 2% effect not determined by 
the theory). At large positive a-r neglecting the exponentially small contributions to U, 
one obtains: 
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a T 



8 2y/2 



— + 2 + 2arc sinh 
8 

16 



l-21og2 + 21oga T + ^- -j- + 



(It 
.TJ 

320 

tjCLrp 



(45) 



On the other hand, the high temperature expansion of the optimized gaussian is 

9 4 18 1324 . , 

21og47r 2 + 21oga T + — - — + 7T7?r (46) 



Cbrp Uirp C/Cbrp 

One observes that the high temperature expansion of two theories are in remarkable 
agreement up to the order -^. 



C. Magnetization, 2D 

Experiments on great variety of layered high Tc cuprates (Bi or Tl || based) show that in 
2D, magnetization curves for different applied field intersect at a single point (M*,T*). The 
range of magnetic fields is surprisingly large (from several hundred Oe to several Tesla). 
Assuming this it is easy to derive the scaled LLL magnetization just from the existence of 
the point. The dimensionless LLL magnetization is defined as p6 | 

dfeffM 



miOTj 



da 



T 



and the measure magnetization is 

e * h 
M = 



m ab 



<H 2 > 



e*h 
m ab 



\A 




(47) 



(48) 



where ip is the order parameter of the original model, and ip r is the rescaled one, which is 
equal to %M Thus 



M 



e*h buj 



-mlaT)- 



m ab V 47r 
Using the definition of a? = —r] 1 ~t- b ,ri = (27t 2 Gi) ,b can be written as 
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(49) 
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Thus eq.( 



implies that 
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If we assume that the experimental observation that all the magnetization curves intersect 
at some point (T*, M*), m (ar) is 



m (o T ) =C 1 [a T ± JC 2 + 4 ; (52) 



m ab M* „ t 2 l-t* 



Ci= " ^ C 2 = Arr 



e*h\l-t*\ t* 

On the other hand, if we require that the first two terms of the high temperature 
expansion of eq.([52"|) and the high temperature expansion of the magnetization are equal, 
one finds that. 

C l = l - ;C 2 = 16 

When we plot this line on Fig. 7 (the dotted line) we find that at lower temperatures the 
magnetization is overestimated. On the other hand magnetization of the theory of 
Tesanovic et al (the dashed line on Fig. 7) underestimate the magnetization. The OPE 
results are consistent with the data within the precision range till the radius of 
convergence a>T = —3. It is important to note that deviations of both the phenomenological 
formula eq.fl5*2"|) and the Tesanovic's are clearly beyond our precision range. 
We conclude therefore that although the theory of Tesanovic et al is very good at high 
temperatures (deviations only at the order 1/af*) they become of the order 5-10% at 
a>T = —3. The advantage of this theory is however that it interpolated smoothly to the 
solid and never deviates more than 10%. The coincidence of the intersection of all the lines 
at the same point (T*, M*) cannot be exact. Like in 3D it is just approximate, although 
the approximation is quite good especially at high magnetic fields. 



D. Specific heat, 2D 

Specific heat OPE result is compared on Fig. 8 with Monte Carlo simulation of the same 
model by Kato and Nagaosa JT5J (black circles), the phenomenological formula following 



from eq.(^) (dotted line) and the theory of Tesanovic et al |n| (the dashed line). The 
agreement with the direct MC simulation is very good. 

E. Magnetization in 3D 

We compare here our results on the LLL scaled magnetization with the Monte Carlo 
simulation of the LLL system by Sasik and Stroud |T5| . They are actually more precise in 
3D. The Fig. 9 contains several OPE approximants (n = 0, 1, 2, 3, 4, 7, 8) and their data on 
all three magnetic fields (representing 2T, 3T and 5T in model YBCO). According to the 
criterion of the "best root" the best approximant should be n = 7. Clearly up to the radius 
of convergence the agreement is within the expected precision. 
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V. CONCLUSION 

In this paper we obtained the optimized perturbation theory results for the LLL Ginzburg 
- Landau model. 

to both the 2D and 3D LLL model. It allows to obtain a convergent series (rather than 
asymptotic). The magnetization and specific heat of vortex liquids with definite precision 
are calculate. One the basis of this one can make several definitive qualitative conclusions. 
The intersection of the magnetization lines in only approximate not only in 3D (the result 
already observed in Monte Carlo simulation jllj), but also in 2D. The theory by Tesanovic 



|T3T| , which uses completely different ideas, describes the physics remarkably well in high 



temperatures and deviates on the 5-10% precision level at ax = —2 in 2D . 
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Figure captions 



Fig. 1 

Feynman rules. Fig. la, Fig. lb, Fig. lc are propagator, four line vertex, mass insertion 
vertex. 

Pig. 2 

Feynmann diagrams for f n [G] for n = 0. 

Fig. 3 

Fig. 3 and Fig. 2 are Feynmann diagrams for f n [G] for n = 1. 

Fig. 4 

Fig4, Fig. 3 and Fig. 2 are Feynmann diagrams for f n [G] for n = 2. 

Fig. 5 

Summing all bubble diagrams. 

Fig. 6 

2D OPT energy and non-optimized energy at different orders (denoted by numbers and T 
plus numbers respectively). One can see clearly OPT series are convergent, for example at 

fig. 7 

2D Magnetization plot: experimental data (Jin. et.alm ref. ||), magnetization from OPT 
for n = 10 (line) , from phenomenological formula (dotted line) and Tesanovic' theory 
(dash-dotted line) 

Fig. 8 

2D specific heat plot. Monte Carlo data by Kato and Nagaosa in ref. ||15|| , specific heat 
from OPT for n = 10 (line) , from phenomenological formula (dotted line) and Tesanovic' 
theory (dash-dotted line) 

Fig. 9 

3D Magnetization plot. Monte Carlo data by Sasik and Stroud in ref. [ 15j , specific heat 



from OPT of different orders denoted by numbers, n = 4, 7 are best approximants among 
them. 
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